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Abstract 

Sequential dynamical systems (SDS) are used to model a wide range 
of processes occurring on graphs or networks. The dynamics of such 
discrete dynamical systems is completely encoded by their phase space, 
a directed graph whose vertices and edges represent all possible system 
configurations and transitions between configurations respectively. 

Direct calculation of the phase space is in most cases a computa- 
tionally demanding task. However, for some classes of SDS one can 
extract information on the connected component structure of phase 
space from the constituent elements of the SDS, such as its base graph 
and vertex functions. 

We present a number of novel results about the connected com- 
ponent structure of the phase space for /c-threshold dynamical system 
with binary state spaces. We establish relations between the structure 
of the components, the threshold value, and the update sequence. Also 
fixed-point reachability from garden of eden configurations is investi- 
gated and upper bounds for the length of paths in the phase space are 
shown to only depend on the size of the vertex set of the base graph. 

1 Introduction 

Graph dynamical systems (GDS) constitute a natural framework to model 
distributed systems such as biological networks or epidemics over social net- 
works. Common to all types of GDS are their defining elements: a finite 
graph Y with vertex set v\Y] = 1, 2, . . . , n and a finite set K whose elements 
are associated to each vertex v of Y. In addition, the transition of each 
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vertex state at each time step depends on the states of its nearest neighbors 
through a vertex function. The actions of the single vertex functions are 
cahed vertex updates. A system update over the entire vertex set is deter- 
mined by a scheme specifying a mechanism by which the individual vertex 
functions are sequentially applied so as to induce a dynamical system with 
map F : ^ K". 

Originating from applications of computer science and theory of simu- 
lations, [2, 3, 4] sequential dynamical systems (SDS) generalize the concept 
of cellular automata (CA) [6, 5] by employing arbitrary finite graphs with 
no restrictrions on the vertex functions (in contrast, CA are usually defined 
on lattices with one single vertex function) and most importantly by being 
characterized by an asynchronous update mechanism specified, for instance, 
by an ordering of the vertex set. 

Applications of SDS include concurrent processes, distributed protocols 
[12], approximate discrete event simulations [7], transport simulations on 
irregular grids [9, 10], image treatment and reconstruction [8], and the mod- 
eling of the mechanisms of gene transcription in messenger RNA [11]. 

The phase space associated to a SDS map F : — >■ K"" is the finite 
directed graph with vertex set K"^ and directed edges {x,F{x)). The struc- 
ture of the phase space is governed by the properties of the base graph Y, 
the vertex functions /„, and the update scheme. The phase space encodes 
the entire dynamics of the system by representing all system states and all 
possible transitions among them. The fundamental properties of SDS have 
been studied in [20, 21, 22, 23, 24, 12]. 

Understanding the long term dynamics of a SDS and being able to com- 
pare the evolution of different system configurations are vitally important 
tasks in many situations. However, the exhaustive calculation of the phase 
space is often resource demanding. Indeed, even for a sequential dynamical 
system over a relatively small base graph, the brute force calculation of peri- 
odic points and transient configurations is a computationally hard problem 
[12, 15, 1, 18, 19]. Research in SDS therefore focuses on deducing the prop- 
erties of the phase space from the features of the constituents of the system. 
This analysis has a local to global character. 

The purpose of the present article is to investigate the connected compo- 
nent structure (therefore unravelling the dynamics of the system completely) 
of fc-threshold SDS with binary states over different classes of base graph. 
Threshold SDS arise naturally in the context of gene annotation methods 
such as functional linkage networks [16]. Furthermore, this class of graph dy- 
namical systems is particularly well-suited for reachability studies [17] since 
they exclusively admit fixed points as periodic points. 

We present a collection of results concerning the dependence of the com- 
ponent structure of different phase spaces on the value of the threshold and 
base graph. The base graphs we considered are complete graphs, circle 
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graphs, and two simple examples of tree: stars, and line graphs. A case 
of particular interest occurs when the threshold value is set to 2. In this 
scenario phase spaces consist of star-shaped components centered at a fixed 
point, and isolated fixed points, and this feature is shared by all four types 
of base graph. This special situation greatly facilitates the determination of 
basins of attractions and w-limit sets. 

In addition, we find upper bounds for the length of directed paths within 
connected components and, construct update orders for which paths of max- 
imal length occurr. Our upper bounds refine the results presented in [15] for 
the cases in question. 



2 Preliminaries 

Let y be a simple, undirected graph (henceforth just a graph). Let v[Y] = 
{1,2,..., n} be the vertex set of Y and e[Y] be its edge set. 

Given v € v[Y] we use n[v] to denote the set of vertices that are contained 
in the 1-neighborhood of v ordered increasingly. If d{v) indicates the degree 
of V, then \n[v] \ = d{v) + 1. 

Let K he a finite set of states. To each vertex in v[Y] we associate a 
vertex state G K. The n-tuple y = {yi,y2, ■ ■ ■ , y-n) is called a system state 
or system configuration. The restriction of y to the vertices of the set n[v] 
is referred to as the state neighborhood of v and is written as y[v]. 

The state of each vertex is updated through a vertex function fy : 
K defined on y{v\. Vertex functions can be incorporated into Y -local Junc- 
tions Fy : K"" ^ K"- as 



FY{yi,y2,---,yn) = {yi,y2, ■ ■ ■ ,yv-i, fv{n[v]),y^+i, . . . ,yn). (l) 
We denote the family of Y-local functions over the vertices of Y by Fy = 

{Fv)vev[Y]- 

Let Sy be the set of permutations (the symmetric group) over v\Y]. 

Definition 2.1. Let y be a graph, K he a finite set, vr = (vj^ ,Vi^, . . . ,Vi^) G 
Sy he a permutation over v\Y], and Fy the family of Y-local functions 
associated to the vertices of Y. We say that the triple (y, Fy,7r) is the 
permutation sequential dynamical system (which we simply refer to as SDS) 
over the base graph Y with state space K and update sequence tt. The 
associated dynamical system map (SDS map) is Ft^ = K^, 



F^ = F,,„oF,,_^o...oF„,^. (2) 

Consider a sequence of maps {gi)f^i with gi : ^ K. We say that 
the maps gi induce the family of vertex functions (/i;)t,g«[y] by setting = 
9d{v)+i- The SDS obtained from these vertex functions is called the SDS 
over y induced by the sequence {gi)f^i- 
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Accordingly, in the case of regular graphs, an SDS is induced if all vertices 
of Y have identical vertex functions induced by gi. 

The phase space of an SDS map is the directed graph r(F7r) with 
vertex set and edge set 

e[r(F^)] = {(y,F.(y)) lyeK"}- 
Starting from a system state x, one can construct the following sequence 

0+(x) = (x,F.(x),F2(x),F3(x),...) = (F^(x))r=o- 

The sequence 0~^{x) is called the forward orbit (or a time series) of 
x under Ft^. Since we only consider finite base graphs and vertex states 
are taken from a finite set, all forward orbits contain a finite number of 
distinct vertices. By definition, the orbit of a given system state traces out 
a directed walk through the phase space. The direction is from the initial 
vertex x onward, that is, we have the alternating sequence of vertices and 
edges 

(x, xF,(x), F„(x), F^(x)f2 (x), f2 (x), . . . ). 

Each of such walks is connected since it contains a path between any two 
of its vertices. Since any vertex of r(F7r) can be used as the initial point 
of an orbit, the phase space of an SDS map is the collection of all forward 
orbits, which its vertices generate. 

Given a system state y the predecessors of y are those configurations x 
such that F7r(x) = y. A configuration that does not admit predecessors is 
called a garden of Eden (GOE). Such configurations can only occur as initial 
states of an SDS and can never be generated during its evolution. 

As in the continuous case, one can recover the notions of periodic points 
and attractors for graph dynamical dynamical systems. 

A fixed point of an SDS is a system state y for which F,r(y) = y- We 
say that an SDS cycles through a (finite) sequence of states (yi,y2, • • • ,yr) 
if F^(yi) = y2, F^{y2) = yg, . . . , F^)(y.J = y^. The system states in 
(yi,y2, . . . ,yr) are called periodic states. Thus, a fixed point is a periodic 
cycle with one single system state. 

Due to the fact that for finite base graphs (the only ones that we will 
consider in our study) the phase space of an SDS map possesses only a finite 
number of nodes and that each node of the phase space has out-degree 1 (SDS 
are deterministic) , the forward orbit of any system state will eventually reach 
a periodic cycle [14]. 

Given an SDS map F^ : ^ K", the set M C for which F^(M) C 
M is called an invariant set of F^r. An tj-limit set of a state of is the 
set of periodic points of F^r, z € i^" such that F™(y) = z for some m > 0. 
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The basin of attraction of an invariant set M C K"^ is the set of the points 
of K"' which possess w-hmit sets contained in M. 
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Figure 1: a) The graph Circ4. The phase spaces of the SDS generated by 
the vertex functions norf:{xi,X2, . . . , x^) = (l + xi) x (1 + X2) x • • • x (1 + 
over Circ4 for the update sequence (0, 2, 1, 3) (b)) and (0, 1, 2, 3) (c)). 

We set sum„ : Jf" N° to be 



sumg(yi, y2, • • • , Vq) = ^yi (computed in N° 



i=l 



Definition 2.2. For any integer k the k-simple threshold function t^ : 
F2 is defined by 



if sumq(y) < k 

1 otherwise 



We wiU use the term k-threshold SDS to refer to an SDS induced by 
A;-simple threshold functions. 

We exclusively consider /c-threshold SDS with vertex states in the Boolean 
set F2. To avoid repeated mentioning of their various constituents we will 
refer to this special type of /c-threshold SDS as binary k-threshold SDS. 

Any threshold SDS only admits fixed points as periodic points [12, 14]. 
This characteristic feature makes this class of SDS particular well-suited for 
reachability [17] and perturbation [14] studies. 

The main goal of our analysis is to investigate the structure of the 
phase space of binary k-threshold SDS and to understand how the prop- 
erties and size of their connected components vary as we change the base 
graph, the value of the threshold, and the update sequence. An estimation 
of the "length" of a connected components can be carried out through the 
fixed point reachability of its GOE configurations, that is, the length of the 
GOE-to-fixed point paths. In other words, given the SDS (y, T'^,7r) we will 
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evaluate upper bounds and exact values of the smallest integer r such that 
(T^)'-(y) = (T^)^-+i(y) for the GOE y [15]. 

We use the standard notation introduced in [12, 13] for the definition 
and vertex enumeration of complete graphs Kn, star graphs Star„_i , circle 
graphs Circ^, and line graphs Line„ over n vertices. 

3 Phase Space Structures 

In this section we investigate how the topology of the connected compo- 
nents of the phase space depends on the threshold value. The system states 
(0,0,..., 0) and (1,1,..., 1) will play a special role in our investigation, and 
for sake of simplicity we will denote them by and 1 respectively. 

The trivial threshold A; = clearly makes little sense in the context of a 
binary state space (every possible system state is updated to 1 after just one 
application of the SDS-map for any base graph and any update sequence), 
and we will discard it altogether in our analysis. 

In addition, we always assume that 1 < A; < d{v) + 1 where d{v) is the 
greatest degree of a vertex in Y. In fact, if A; > d{v) + 1 no vertex state can 
ever be updated to 1 and any system state is updated to for any update 
order. 

Under these assumptions, binary /c-threshold SDS always admit at least 
one fixed point. 

Lemma 3.1. For a binary k-threshold SDS, (K, (r^)^,7r) with k > 1, the 
state is a fixed point of the SDS-map T^. 

Lemma 3.2. For a binary k-threshold SDS, {Y, (T^)t,,7r) the system state 1 
is a fixed point of the SDS-map T^, if for all vertices v € v[Y], d{v) + 1 > A;. 

We will treat the case of 2-threshold SDS separately. Not only is 2 
the first nontrivial threshold value for the binary state space, it is also the 
value for which the phase spaces over the four types of base graph possess 
connected components with the same structure. 

For the case of 1-threshold SDS we make the following observation. 

Proposition 1. After exactly one system update all the vertex states in a 
neighbourhood containing at least one 1 state when the update is applied, will 
be changed to 1. This happens independently of the update sequence. 

3.1 Complete Graphs 

Complete graphs Kn offer the advantage of each vertex possessing a "global 
neighborhood", that is, it has all the others vertices in v[i^n] as nearest 
neighbors. This feature simplifies considerably the analysis of the phase 
space of threshold SDS over complete graphs. 
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Theorem 3.3. Let [Kn, (T^^)^,7r] be a binary k-threshold SDS with 1 < k < 
n. Then the phase space of 1^ consists of the union of 2 star-shaped con- 
nected components, whose centers are the fixed points and 1 . In addition, 
the structure of the phase space is independent of the update sequence vr. 

Proof. Since every vertex of Kn belongs to a global neighborhood, the update 
of a vertex state affects the update of all the other vertex states. For a 
threshold function the update of Xi depends on sum„(x[z]), so clearly the 
order in which the vertices are updated is unimportant to the outcome of 
the sum. 

Given an initial system state x € F2, the system will evolve in either of 
the following ways. 

• If sum„(x) < k: the initial state contains less than k vertex 1 states, 
and each of the vertex states will be updated to 0. The update ends 
with the fixed point 0. 

• If sum„(x) > k: the initial state contains more than k 1 states and 
all vertex states will be updated to 1. The update ends with the fixed 
point 1. 

So the orbit of any system state is attracted either by or 1 after one 
single iteration of the SDS-map. 

□ 

Counting GOE configurations and determining the basins of attractions 
of and 1 thereby become straightforward. 

Corollary 1. Let [Kn, (r^')„, vr] he a binary k-threshold SDS with 1 < k < n. 
The basin of attraction of 0, B{0), consists ofJ2i=i CI) points. 

Proof. We have (") system states that possess i vertex 1 states. To obtain 
the number of system states that end up at after one update, we sum over 
i from 1 to A; — 1 to cover all the possible system states with the number of 
vertex 1 states ranging from 1 to /c — 1, that is, |5S(0)| = X)i=i O- 

□ 

Corollary 2. Let [Kn, (r^)«, vr] be a binary k-threshold SDS with 1 < k < n. 
The basin of attraction B{1) of 1 consists of 2" — J2^=k il) ~ ^ points. 

Proof. All the system states x ^ B{0) that are not or 1 belong to 
Thus, we subtract |iB(0)| = X]i=i (1) from the total number of possible 
system states, 2", then we subtract 2 to account for the two fixed points, 
and 1. We thus obtain |^(1)| = 2" - J2i=i Q "2. □ 
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3.2 Star„ 

The center of a star graph occupies a special position, not only is it the 
only vertex with a global neighborhood, it is the only vertex common to all 
neighborhoods of the vertices in Star^. With reference to the shape of a star, 
we introduce the following short-hand terminology. We refer to the elements 
of f [Star„] that are not the center as "arm vertices" and denote their states 
as . The state of the center vertex is then called Xr- 
For 1-threshold system one has 

Theorem 3.4. Let [Si:ar„, (T^)„, vr] be a binary 1-threshold SDS. Then the 
phase space of consists of a directed tree rooted at the fixed point 1 for 
which the leaves have a distance at most 2 from the root, and the isolated 
fixed point 0. 

Proof. We observe that all the vertex states that are 1 or have neighboring 
vertices with state 1 will either stay 1 or be changed to 1. We have three cases 
depending on the state of the center and the 1 states in its neighborhood. 

• Xc = and sum„+i(x[c]) = when Xc is updated. 
This situation corresponds to the fixed point 0. 

• Xc = and sum„+i(x[c]) > 1 when Xc is updated. 

The states Xa that are updated before Xc maintain their states since Xa 
is the only contribution to the sums sum2(x[a]). Then Xc is updated 
to 1 and so is every vertex state updated after Xc- Since now Xc = 1 
all the states Xq that are not 1 already will get updated to 1 at the 
next system update. The system state will be attracted by the fixed 
point 1 after two systems updates at most. 

• Xc = 1 when Xc is updated. 

Every vertex state gets updated to 1 independently of the update se- 
quence since Xc always equals 1. We reach the fixed point 1 after 
exactly one system update. 

□ 

The following corollary is a consequence of the proof of theorem 3.4. 

Corollary 3. The SDS-map possesses only two fixed points, the isolated 
fixed point and the fixed point 1. 

Any system state of a binary 2-threshold SDS over Star„ is fully charac- 
terized by the position of the center vertex c in the update sequence vr and 
by its state Xc- This observation allows us to devise four simple and very 
useful rules for 2-threshold SDS over star graphs. 
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Proposition 2. Any state Xa will be updated to when Xc = 

Proposition 3. If Xc = 1, the states Xa maintain their values under vertex 
updates. 

Proposition 4. To update Xc from to 1, we must have at least two 1 states 
in Xa in x[c\ when Xc is updated. 

Proposition 5. To keep Xc = 1 upon update, x[c] must contain at least one 
state Xa = 1. 

Theorem 3.5. Let [Starn, {T^)v,'^] be a binary 2-threshold SDS. Then the 
phase space of consists of star-shaped connected components centered at 
fixed points. 

Proof. The rules above allow us to fully characterize the phase space of T^ 
through the following five cases: 

• Xc = and there are less than two states = 1 when Xc is updated. 

If all states Xa are 0, by proposition 2, we have the fixed point 0. 
Otherwise there is one single = 1. Independently of the sequence vr, 
Xc can never be updated to 1 as only one vertex would contribute to 
sum„+i. In addition, by proposition 2 the arm vertex with Xa = 1 will 
be updated to Xa = 0, while all the other vertices retain their state. 
Since Xc = before it is updated, all the arm vertices that were updated 
before c were necessarily updated to according to proposition 2. 

After one system update we therefore recover the fixed point 0. 

• Xc = and there are two or more states Xa = 1 when Xc is updated. 

All the vertices that are updated before c will get the state by propo- 
sition 2. Since sum„_|_i(x[c]) > 2, Xc gets updated to 1 (proposition 4). 
The states that were equal to 1 when Xc was updated will be left 
unchanged when they are updated (proposition 3). The next system 
update leaves the system state unchanged since Xc = 1 ensures that 
all the vertices updated before c maintain their state (propostion 3). 
Again x[c] contains more than two states equal to 1, leaving Xc = 1 
(proposition 4), so that the states that are updated later also will not 
change. 

• Xc = 1 and all the x^ = when Xc is updated. 

If all Xa = 0, they keep their states (proposition 3) and Xc gets up- 
dated to too, independently of the update sequence since sum,i+i(x[c]) = 
1. Again, we recover the fixed point 0. 
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• Xc = 1 and there is one Xq = 1 in x[c] when c is updated. 

All arm states updated before c retain their initial value by proposition 
3. The center state Xc stays 1 since sum„_|_i(x[c]) = 2. By proposition 
3 all the arm states that are updated after c also maintain their values. 
Such configurations are therefore fixed points. 

• Xc = 1 and there are two or more than two 1 states in x[c] when c is 
updated. 

All the vertices updated before c maintain their original state and so 
does c when it is updated (propositions 3 and 5). The state Xc = 1 en- 
sures that all the vertices updated after c keep their state (proposition 
3). Such system state configurations are therefore fixed points. 

We have thus proven that in each of these cases we have (T^)^(x) = 
T^(x) for all system states x € F2+^ □ 



00101 00100 00011 00010 00001 00111 01011 01101 01110 




Figure 2: a) The graph Star4. b) Phase space of the 3-threshold SDS over 
Star4 with the identity update sequence. 

Corollary 4. The number of fixed points in the phase space of is |F/x(r^)| = 
2". 

Proof. From the proof of theorem 3.5 we know that, apart from 0, the fixed 
points of are characterized by Xc = 1 and by sum„_|_i(x[c]) > 2. There are 
n ways to construct a fix point with only one x^ = 1. Requiring two states 
Xa to be equal to 1 allows for the construction of (2) fixed points. Therefore 
there are {^) fixed points containing k vertex states Xa = 1. Counting 0, we 
have 

1=1 ^ / 1=0 ^ 

□ 
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Prom the proof theorem 3.5, we see that n fixed points correspond to 
the configurations with one Xq = 1 and Xc = 1. We observe that the results 
of theorem 3.5 and corollary 4 hold for 2-threshold SDS only. For a generic 
fe-threshold system with k > 2, the phase space consists of a single connected 
component, a tree rooted at the fixed point 0. 

Theorem 3.6. Let [Star„, (T^')^,, vr] be a binary k-threshold SDS with 2 < 
k < n + 1. Then the phase space of consists of one single connected 
component, a tree rooted at for which the leaves have a maximum distance 
of 2 from the root. 

Proof. Notice again that the state Xc and the position of the center vertex 
c in the update sequence vr fully characterizes the phase space ofT|^. None 
of the vertices on the arms of the stars can ever be updated to = 1 since 
sum2(x[o]) < 2 always, while we need to have at least k vertex 1 states in 
x[c] to update Xc to 1. We basically have two scenarios: 

• sumfi+i(x[c]) < k when c is updated. 

Then Xc is set to along with all the other vertex states of the star 
graph. The update creates the fixed point 0. Such points have unit 
distance from 0. 

• sum„+i(x[c]) > k when c is updated. 

During the system update, the state Xc is set to 1 while all the other 
vertex states become 0. The next update will set the system state to 
0. The system passes through the state in which only Xc = 1 before 
reaching upon a second update. Such points have distance 2 from 0. 

We can thus conclude that the phase space of is constituted by a single 
tree-like (i.e. a directed tree with a loop at the root) connected component 
rooted at 0. The leaves have distance either 1 or 2 from the root. 

□ 

From the proof above it follows immediately that 
Corollary 5. The system state is the only fixed point for the map T^. 

3.3 Circ„ 

Consider now 1-threshold SDS over circle graphs (or n-gons). We notice 
that (by proposition 1) all system state configurations that contain at least 
one vertex 1 state will eventually be updated to 1 upon iterative application 
of the SDS-map. We can thus conclude that all the system states but 
belong to the same connected component rooted at the fixed point 1. As 
a consequence of proposition 1 for a vertex 1 state to be created after a 
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system update, at least one of its neighboring states must necessarily be 1 
before the SDS-map is applied. System states that possess configurations 
consisting either of one single vertex 1 state or where all the 1 states are 
separated by one or more state are therefore GOE configurations. 

Lemma 3.7. The phase space of a binary 1-threshold SDS over CirCn is 
constituted by the isolated fixed point and a tree-like connected component 
rooted 1. 

With this in mind, we construct an update sequence for which the GOE- 
to-root path in the component containing 1 is the longest possible. We 
begin by considering a system state for which only the i-th vertex state is 
different from 0. Such configurations are the GOE with the greatest number 
of vertex states that can be updated to 1. At least two vertex states must 
be updated to 1 at each application of the SDS-map (if only one state were 
updated to 1, we would have that a neighboring a 1 is unaffected by the 
system update, thus contradicting proposition 1). The update order that 
minimizes the change in vertex states at each application of the SDS-map 
(thus achieving the longest GOE-to-root path) is necessarily one for which 
two states are updated to 1 at each application of the SDS-map. 

Starting from the vertex farthest away from i we update the vertex states 
on the sides of (i + \ n/2\ ) mod n alternatively. We have the update sequence 

vr ={{i + \n/2\) modn, [i + \n/2\ — 1) modn, {i + [n/2j -|- 1) modn, 
{i + \n/2\ — 1) mod n, {i + [n/2j + 2) mod n, . . . ,i). 

As the SDS map is applied iteratively more and more 1 states are added 
to the sides of Xi until the system reaches 1 after [n/2\ updates. 

Theorem 3.8. Let (C/rCn, (TJ)^,, vr) be a binary 1-threshold SDS. There ex- 
ists at least one update order fr for which the phase space of T| contains 
a path of length [n/2j and this is the maximum GOE-to-root path length 
attainable in the phase space. 




Figure 3: Compatible neighborhoods on Cirq 



We continue by considering 2-threshold SDS over CirC; 
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Theorem 3.9. Let [CirCn, {T^)v,tt] be a binary 2-threshold SDS. Then the 
phase space of consists of star-shaped connected components centered at 
fixed points. 

Proof. By construction of Circ„, for each vertex i the neighborhood state 
x[i] = Xj, Xj+i) has two overlapping neighborhoods x[i — l] = (xj_2, Xj) 

and x[i + l] = (xj, Xj+i, with the identifications n + 1 = 1 and n + 2 = 2. 
We will refer to such set of neighborhoods as compatible neighborhoods. 
Compatible neighborhoods are all affected by the update of the vertex state 

We will show that after one system update none of the vertex states in 
x[i] will change again when a second update is run. 
We have two cases: 

• Suppose that the Y-local function T^ is applied to x[i] and x-i is up- 
dated to 1. 

For Xi to be updated to 1, x[i] must have one of the following configu- 
rations at the time of the update: (1, 1, 1), (1,0, 1), (1, 1, 0) or (0, 1, 1). 
After the Y- local function T^ is applied we will have (1,1,1), (1,1,1), 
(1, 1, 0) or (0, 1, 1) respectively. These are all local fixed points of T^. 
In addition, the configurations (1, 1, 1) and (0, 1, 1) "lock" the update 
of Xj+i (if Xj+i were not updated before Xi) to 1 (by compatibility 
there are two 1 state in x[z-f 1]). While the configurations (1, 1, 1) and 
(1, 1,0) lock Xi-i to 1 in the same way. 

In any of these cases there will always be at least two 1 states in x[i] 
when a second update is applied. Thus, the state Xi will remain fixed 
under further system updates and this holds for all the vertices of Circ„. 

• Suppose that the y-local function is applied to x[i] and Xi is up- 
dated to 0. 

For Xi to be updated to 0, x[i] must have one of the following configu- 
rations at the time of the update: (0, 0, 0), (0, 1, 0), (1, 0, 0) or (0, 0, 1). 
After the F-local function T^ is applied we will have (0, 0, 0), (0, 0, 0), 
(1,0,0) or (0,0,1) respectively. These are all local fixed points oiTf. 
In addition, the configurations (0,0,0) and (1,0,0) "lock" the update 
of Xj+i (if Xj+i were not updated before Xi) to (by compatibility 
there are two state in x[i -|- 1]). While the configurations (0, 0, 0) and 
(0, 0, 1) lock Xi^i to in the same way. 

Any of these cases will ensure that there are at least two states in 
x[z] when a second update is applied. Thus, the state Xi will remain 
fixed under further system updates and this holds for all the vertices 
of Circ„. 
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3 PHASE SPACE STRUCTURES 



Any given system state is therefore either a fixed point or wih reach a 
fixed point after exactly one update. 

□ 
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Figure 4: a) The graph Circs. b) Phase space of 2-threshold SDS with 
5-gonal base graph and the identity update sequence. 



Although n-gons are the simplest case of circle graphs, there is no im- 
mediate way to extract information about connected components or fixed 
points of a 2-threshold SDS by simple inspection of the system states. The 
most suitable method to obtain all the fixed points of the system is to con- 
struct a local-fixed-point graph as described in [12]. 

When we raise the threshold value to 3, we notice that for a vertex state 
rEj = 1 to remain 1 after the y-local function is applied, all the states in x[i] 
must be 1 as well. Clearly, the fixed point 1 is the only system state for 
which this condition is satisfied for all vertices. If any of the vertex states 
in a system state is 0, the vertex state will propagate within the system 
configuration until is reached after a number of system updates. Graphi- 
cally we have the opposite situation to the case of 1 threshold systems and 
an equivalent version of proposition 1 for which and 1 vertex states are 
exchanged. It follows that for 3-threshold SDS over n-gons GOE configura- 
tions are those containing either one single vertex state or where vertex 
states are separated by one or more 1 states. 

Lemma 3.10. The phase space of a binary 3-threshold SDS over CirCn is 
constituted by the isolated fixed point 1 and a tree-like connected component 
directed towards the root 0. 

Considering the same update sequence constructed for theorem 3.8, with 
the difference that the starting GOE configuration has all vertex 1 states 
except for the one in the i^^ position, we obtain the maximum path lenght 
in a connected component. 

Theorem 3.11. Let {CirCn, (T^)in'^) be a 3-threshold SDS. There exists at 
least one update order tt for which the phase space of contains a path of 
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3.4 Line, 



3 PHASE SPACE STRUCTURES 



length \n/2\, and this is the maximum path length attainable in the phase 
space. 

3.4 Line„ 

If we set the threshold value to 1 for threshold SDS over line graphs, one 
recovers the same situation as for 1-threshold SDS over n-gons with the 
obvious difference in the neighborhoods of the end vertices. Again, we have 
that it is sufficient to have a single vertex 1 state to ensure convergence of 
the system state to 1 upon iterative update of the system for any update 
sequence as proposition 1 holds for any 1-threshold SDS. Analogously, GOE 
configurations will be those system state containing one single vertex 1 state 
or more vertex 1 states separated by at least one state. 

Theorem 3.12. The phase space of a binary 1-threshold SDS over a line 
graph is constituted by a the isolated fixed point and a tree-like connected 
component rooted at the fixed point 1 . 

We can therefore apply the same arguments we have used for n-gons 
to find an update sequence that gives the longest GOE-to-root path in the 
connected component containing 1. 

Starting at the GOE state (0, . . . , 0, 1) we use the identity update se- 
quence. It will take exactly n — 1 updates to reach the 1. This is the 
maximum number of applications of the SDS-map, since to produce at least 
n steps one should have an update order that does leave the entire system 
state unchanged under at least one update. Clearly this is only possible for 
fixed points. 

Theorem 3.13. There exists at least one update sequence that maximizes 
the GOE-to-root path length in the phase space of a binary 1-threshold SDS. 
The maximum length attainable is n — 1. 

Besides being important to numerous applications, line graphs constitute 
an example of how our results can be combined to work out the phase space 
component structure of new SDS from old ones. Theorem 3.5 and theorem 
3.9 allow us to derive the following result. 

Theorem 3.14. Let (Linen, {Ty)v,'^) be a binary 2-threshold SDS. Then the 
phase space of consists of star-shaped connected components centered at 
fixed points. 

Proof. First, one has to distinguish between two type of vertex present in a 
line graph. We have internal vertices, that is, all the vertices v with d{v) = 2 
and the two leaves at the ends of the line with degree 1. 

For the internal vertices we have the same neighborhood structure that 
we have seen in the proof of theorem 3.9. Thus, the result for 2-threshold 
SDS over Circ„ apply for internal vertices. 
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4 CONCLUSIONS 



The update of the leaves depends only on the state of the internal vertex 
to which they are adjacent, in other words it behaves exactly as the update 
of an arm vertex in a star. Thus, theorem 3.5 ensures that the vertex states 
of the extremities are fixed after exactly one system update. 

So for any vertex v G Line^, x['y] reaches a "local fixed point" after at 
most one application of T^. 

We conclude that for any system state (T^)^(x) = T^(x) holds. □ 

Again, fixed points of line graphs can be determined through a general- 
ization of the local fixed point method in [12]. The basins of attraction can 
be determined accordingly. 

For a threshold value A: = 3, we see immediately that the only fixed point 
is 0. In fact, all the systems states containing vertex 1 states will eventually 
be updated to since there is always at least one neighborhood x[v] for 
which sum^xlv] < 3. The GOE configurations for 3-threshold SDS over line 
graph will then have the same characteristics as those for 3-threshold SDS 
over n-gons. 

To set up the update sequence that gives the longest possible GOE-to- 
root path length, we choose the state 1 as starting point. Since for any 
update sequence the end vertex states will necessarily be update to 0, we 
proceed from the ends towards the center vertices in the same fashion as we 
did for n-gons. We thus consider 

7r = ([n/2j, [n/2j - 1, K2j + 1, Ln/2j -2, [n/2\ +2,...,l,n) 

With this update sequence it takes 1 exactly [n/2] applications of the 
SDS-map to be updated to 0. Since for any two vertex states that are 
separated by a string of 1 states on a line graph each system update must 
produce at least two new vertex states. We have that the sequence that 
we constructed gives the longest GOE-to-root path (for the equivalent of 
proposition 1 to hold). 

We collect our observations in the following theorem 

Theorem 3.15. The phase space of a binary 3-threshold SDS over Linen 
is a single tree-like connected component rooted at the only fixed point, 0. 
In addition, there exists at least one update sequence that attains the upper 
hound of \n/2\ to the GOE-to-root path length. 

4 Conclusions 

In our survey of binary /c-threshold SDS we have completely characterized 
the phase space of systems over complete, circle, line, and star graphs. We 
have proved that the reachability of fixed points from GOE configuration 
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Table 1: Maximal GOE-to-fixed-point path length for the four base graphs. 





Kn 


Circ„ 


Line„ 


Star„ 


k = 1 


1 




n — 1 


2 


k = 2 


1 


1 


1 


1 


k = 3 


1 


LtJ 


Lfl 


2 


k>3 


1 


1 


1 


1 



depends on the threshold value in all cases except for complete base graphs. 
We have explicitly constructed update sequences which maximize the GOE- 
to-root distance. This completely solves the problem of the reachability 
of fixed points [15, 1] from any given system state for the SDS in question. 
The work presented in this article suggests that a similar investigation can be 
extended to binary /s-threshold SDS over generic tree graphs and that the 
GOE-to-root maximum length will again depend on the relation between 
the degree of non-leaf vertices and the threshold value. To obtain similar 
results for /c-threshold SDS over regular graph would be also desirable. This 
especially important in the light of the results on the reduction of some 
classes of SDS to SDS induced by threshold and inverted threshold vertex 
functions over regular graphs presented in [15]. 

It would also be of interest to extend the investigation to state spaces 
containing more elements than F2. 
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